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Abstract 

Recently, it was proposed to use measurements of B®(t) — > £)(*)+£)(*)- and B® — > 
£)(*)+£)(*)- decays to measure the CP phase 7. In this paper, we present the ex- 
traction of 7 using this method. We find that 7 is favored to lie in one of the ranges 
[19.4°-80.6°](+0° or 180°), [120°- 147°] (+0° or 180°), or [160°- 174°] (+0° or 180°) 
at 68% confidence level (the (+0° or 180°) represents an additional ambiguity for 
each range). These constraints come principally from the vector- vector final states; 
the vector-pseudoscalar decays improve the results only slightly. Although, with 
present data, the constraints disappear for larger confidence levels, this study does 
demonstrate the feasibility of the method. Strong constraints on 7 can be obtained 
with more data. 
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Recently, two of us (AD, DL) proposed a method for extracting the CP phase 7 
from measurements of B° d (t) -> dW+dW" and -> £>(*)+£)(*)- decays |T]. The 
technique is quite straightforward. Consider the pseudoscalar-pseudoscalar (PP) 
decay B® — > D + D~ . The amplitude for this decay receives several contributions, 
described by tree, exchange, b — > d penguin and color-suppressed electroweak pen- 
guin diagrams [2J: 

A D = ( T + E + p c ) y* Vcd + Pu V* b V ud + (P t + P° w ) V* b V td 

= (T + E + P c -P t - P° w ) V: b V cd + (P u -P t - PfJ V* b V ud 

= Act e iSct +A ut e^e l5ut . (1) 

Here, Act = \(T + E + P c - P t - P^ w )V; b V cd \, A ut = \(P U - P t - P° w )V: b V ud \, P t is the 
b — > d penguin amplitude with an internal 2-quark, and we have explicitly written 
out the strong phases 5 ct and 5 ut , as well as the weak phase 7. The second line is 
obtained by using the unitarity of the CKM matrix, V* b V ud + V* b V cd + V* b V td = 0, 
to eliminate the Vt b Vu term. The amplitude A D for the decay B® — ► D + D~ can be 
obtained from the above by changing the signs of the weak phases. 

There are three observables which can be obtained from a time-dependent mea- 
surement of this decay: B (the branching ratio), a d i r (the direct CP asymmetry) 
and ai nd i r (the indirect CP asymmetry). In terms of the above parameters, these 
can be written 

B = - (\A D \ 2 + \A D \ 2 ) = Al + Al + 2A ct A ut cos 6 cos , 

a dir = - (\A D \ 2 - \A D \ 2 ) = -2A ct A ut sm5smj , (2) 

a mdir = Im(e- 2i ?A D *A D ) 

= -Al sin 2(3 - 2A ct A ut cos 5 sin(2/3 + 7) - A 2 ut sin(2/5 + 2 7 ) , 

where 5 = 5 ut — 5 ct . Here, f3 is the phase of B^-B® mixing, which has been measured 
in the CP asymmetry in B d (t) — > J/ipK s However, these observables still depend 
on four unknown theoretical parameters: the two magnitudes A c t and Aut, one 
relative strong phase 5, and the weak phase 7. We therefore have three observables, 
but four theoretical unknowns. Thus, in order to obtain weak phase information, it 
is necessary to add some theoretical input 0j. 

This input comes from the decay B d — > D+D~, which receives tree, b — > s 
penguin and color-suppressed electroweak penguin contributions 0: 

A D S = {r + p , )v:bVcs + Kv:bVus + {pl + K C w)Vtlv ^ 

= (t' + p' c - p; - p' E c w ) v: b v cs + {P' u - pi - P£,) v* b v us 

« (T' + P> - PI - P' E C W ) V: b V cs = A' ct e^ . (3) 
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^A D *A D )=B 2 -a 2 



dir ^indir 



(The primes on the amplitudes indicate a b — > s transition.) Here, the last line arises 
from the fact that |V^14 a /V^V^| ~ 2%, so that the piece proportional to V* b V us is 
negligible. The measurement of the total rate for B® — > DfD~ therefore yields A' ct . 
We now make the SU(3) flavour assumption that 

sing c A ct = sing c |(r + ^-P/-P^)K^| = 

A* I (r + £ + p c - p t - pc w )v; b v cd \ ' 1 1 

where sin 8 C is the Cabibbo angle. With this assumption, the knowledge of A' ct can 
be used to give us Act- This in turn means that the three observables in B\ — > D + D~ 
depend on only three theoretical unknowns: A u t, 5, and 7. We can therefore now 
solve for 7 (up to discrete ambiguities). 

The explicit solution for 7 is as follows. We introduce a fourth (non-independent) 
observable, a R . 

a R = Re (t 

= A 2 ct cos 2/3 + 2 At At cos 5 cos(2/? + 7) + A 2 ut cos(2/? + 27) . (5) 

One can obtain a R from measurements of B, adir and a inc nr, up to a sign ambiguity. 
One can then easily obtain 

2 a R cos(2/3 + 27) - a indir sin(2/5 + 27) - B 
A « = ^2^1 ' (6) 

Given the knowledge of 2/?, the assumption in Eq. (@J therefore allows us to obtain 
7- 

The leading-order theoretical error in this technique is given simply by the SU(3)- 
breaking ratio of decay constants /d 3 //d PQ- (There are other errors, such as the 
neglect of the E amplitude in At, but these are expected to be smaller.) This ratio 
has been computed quite precisely on the lattice: /d,//d — 1-20 ± 0.06 ± 0.06 [S]. 
With this value, the theoretical error in this method is rather small, so that 7 can 
be extracted from measurements of B®(t) — > D + D~ and B® — > DfD~. 

Unfortunately, at present data on B®(t) — > D + D~ is unavailable, so we cannot 
apply the above method to this decay. However, the vector-vector (VV) decays 
Bd(t) ~* D* + D*~ and B® — > D* + D*~ have been measured. The method can be 
applied in a similar way to these decays, but with some additional complexity as 
described below. The main theoretical error is given by fo* / Jd* ■ 

A modification of this method can be used when the final state is not self- 
conjugate, as is the case for vector-pseudoscalar (VP) final states jHj- Consider the 
decay B® — > D* + D~ . Following Eq. (JTJ), its amplitude can be written 

A(B° d - D*+D~) = At e l&ct + At e*V 5U< . (7) 

(Although we use the same symbols, the amplitudes and strong phases are not the 
same as those for B® — ► D + D~ .) Now consider the decay of a B® meson to the 
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same final state, B® — > D* + D~. The amplitude for this decay is not simply related 
to that for — > D* + D~ since the hadronization is different: in B^ — > D* + D~ , the 
spectator quark is part of the D~ , while in — > D* + D~ it is contained in the _D*+. 
We therefore write 



D* + D~ 



i<5 c * 



(8) 



The measurement of B®(t) — > D* + D~ still yields three observables, -B, a^r and 
dindir, but now they take more complicated forms. The three observables for the 
D* + D~ final state are: 



B 

&indir 



At + A 2 ct + ZAutAct cos(7 + <5) + ^ + A 2 ct + 2At At cos(7 - 5) 
At + At + ZAutAct cos( 7 + 5) - ^ t - Alt - ZAutAct cos( 7 - 5) 



AtAt sin(-2/3 -2 7 -5 + 5-A) + At At sin(-2/3 - 7 - 5 - A) 
+At At sin(-2/3 - 7 + 5 - A) + AtAt sin(-2/3 - A) , 



(9) 



where 5 = 5 Mt - <P, 5 = <5 Mt - <5 ct , and A = 5 ct - 5 ct . 
For the D + D*~ final state, the observables are 



B = 



Q>dir 
Qindir 



Alt + Alt + 2AtAt cos( 7 + 5) + At + At + 2 At At cos( 7 - 5) 
Alt + A% + 2AtAt cos( 7 + 6)- Al t ~ A% - 2 At At cos( 7 - 5) 



AtAt sin(-2/3 - 2 7 + 5- 5 + A) + AtAt sin(-2/3 - 7 + 5 + A) 

+AtAt sin(-2/3 - 7 - 5 + A) + AtAt sin(-2/3 + A) . (10) 



Considering that /3 has been experimentally determined, the 6 observables are 
written in terms of 8 theoretical unknowns: A ut , At, At, At, 7, 5, 5 and A. 
We therefore need 2 assumptions to extract information. These come from using 
the decays B® — > D* + D~ and B^ — > D^D*~ . The measurements of the branching 
ratios for these decays allow us to extract A' ct and At, respectively With motivation 
similar to that for the PP mode above, we assume that 



sin 6L A 



At 



ct 1 

- 1 - , 



sin # c At 



A 



= 1 . 



(11) 



ct 



With these assumptions, the measurements of the branching ratios for B® — > D* s + D~ 
and 5^ — > D+D*~ give us At and At- We now have 6 observables and 6 theoretical 
unknowns, thus we can solve for 7 (up to discrete ambiguities). The main theoretical 
error is the deviation from unity of /d|//d* and /d s //d in the first and second 
assumption above, respectively. 
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In this Letter we extract 7 via this method, using BaBar and Belle data on the 
VV and VP modes [7| We determine the value of 7 to be in one of the ranges 
[19.4°-80.6°](+0° or 180°), [120°- 147°] (+0° or 180°), or [160°- 174°] (+0° or 180°) 
at 68% confidence level (C.L.), where the (+0° or 180°) represents an additional 
ambiguity for each range. The first of these ranges, [19.4° — 80.6°], is the range that 
is favored by other, external information on the Unitarity Triangle (assuming the 
standard model) ^2JEHE3- As we will see, the constraints on 7 come principally 
from the data on VV modes; at present, the VP decays add only a small amount. 
Note also that if we consider a larger C.L. (e.g. 90%), the constraints on 7 disappear. 
Thus, even though we obtain limits on 7 at 68% C.L., our main purpose here is to 
demonstate the feasibility of the method. With more data, the constraints on 7 
will be correspondingly stronger. Eventually this method can be used to obtain a 
precision determination of 7. 

We begin with the analysis of vector-vector decays. In order to extract 7, how- 
ever, we must make additional assumptions. VV final states come in three transver- 
sity states, 0, ||, and _L. The amplitudes and || are CP-even, while _L is CP-odd. 
The data shows that the D*D* final state is almost entirely CP-even, i.e. the _L 
amplitude is negligible ^T]. Unfortunately, at present experiments cannot distin- 
guish between the and || amplitudes. This has the following effect. For a single 
transversity a, the three observables in B^(t) — ► D* + D*~ can be written similarly 
to Eq. ©: 

B a = A\ a + A\ a + 2 Act, a Aut,a cos 5 a cos 7 , 

a dir = -2 A ct: a A uti(T sm 5 a sin j , (12) 
<ndir = ~At,a sin 2/5-2 Act,, Aut,a cos 5 a sin(2/3 + 7) - A 2 ut ^ sin(2/3 + 2 7 ) , 

where S a = — 5^. However, for all three observables, what is measured is the 
sum of the helicities and ||: 

B = [Alt,o + ^-lt,a + 2 Act^Aut,*, cos <S CT COS7 , 

<T=0,|| 

a-dir = [-2A cti(T A utta sm 5 a sm'j} , (13) 

O-=0,|| 

a-indir = [-"4ct,<7 sin 2/3-2 Act,*, A ut ,a cos 5 a sm(2(3 + 7) - A\ a sin(2/5 + 27) . 

(7=0,|| 

Because the parameters A ctl A ut and 5 are different for the two helicities and ||, 
there are too many theoretical unknowns to apply the method. 

We would like to cast the expressions for the observables in Eq. (|13|) in the same 
form as those in Eq. @. In order to do this, we must relate the parameters for 
the and || helicities. Specifically, we assume that the strong phases are equal: 
5q = S\\ = 5. We also assume that the amplitudes are proportional to one another 
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(with the same proportionality constant, c): Act,o 
our assumptions are: 

, o = 



A, 



ct. 



A 



cAct,\\, Autfl = cA ut ,\\. That is, 
II ■ (14) 



uf, 



With these definitions, the observables in Eq. (j!3|) take exactly the same form as 
Eq. (J2J), and the method for extracting 7 can be applied. 

In fact, these assumptions are theoretically reasonable. The amplititudes for 
B° d D*D* and B° d -> D*D* are given by 



Avv,a = Act <(T e lS « + A u t a e 



A V v,a = A ct „e iS ° + A ut „e l ^e 



(15) 



with a — 0, 
amplitudes: 



being the three transversity states. We now define the ratios of 



k 



ct 



A c t,0 i(S^-Sf) , A u tfl itS^-Sff*) 

= e 11 , k ut = — — -e 11 . 



ct, 



A 



(16) 



Our assumptions are equivalent to k ct = k ut , i.e. both amplitudes and phases are 
equal. Below, we investigate the extent to which these relations hold true. 

We note that, in general, we can write any amplitude in terms of factorizable 
and nonfactorizable pieces: 



ct.F 



Aut,cr A u t a & a 1 "I - s a c 



(17) 



where we denote the factorizable contributions by an index l F\ The quantities r CT , 
p a , s a and A CT parametrize the ratios of nonfactorizable and factorizable amplitudes. 

Consider first only the factorizable contributions. In this case, the ratios of 
factorizable amplitudes are 



h F _ -^ctfl i(6 -5 ) 
^ct — AF 1 

ci >ll 



^4 j. n ■fS'U-t.F cUt,F\ 

h F _ ^utfl j{8 Q ) 
^ut — AF 
ut,\\ 



(18) 



The factorizable amplitude for the decay B d — > D*D* is given by |T2| 



Gi 



A[B° d -> D*D*] a = -j=XP% 



V2' 



(19) 



where X = X\ + X 2 , and 

X 1 = 
x 2 = 

PL* = 



VcbV* d (I2 + a\ + a\ — a\ — a\§ 
VubV ud a>4 + <4 — 04 — , 
m D *fD*£*D* (D*\ ct m (1 - 75)6 



(20) 
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In the above, a,j = Cj + Cj-\/N c , where the Cj are Wilson coefficients. From Eqs. (fTT)|) 
and (J2"U|) . we can read off the individual factorizable amplitudes: 



iS' fT 



Gf 

71 



A* e 1 ' 



Gj_ 

V2 



XoPl 



2 r D« 



(21) 



However, note that the strong phases come from a^1 , and appear only in the factors 
Xi t 2 in Eq. ([20)1 . These factors are independent of transversity. That is, the relative 
strong phases between the factorizable ct and ut amplitudes are independent of the 
polarization state, leading to Arg(A^) = Arg(A^). Furthermore, the expressions in 
Eq. (}2Tj) above lead to 



u t 



Pi 



(22) 



Thus. 



we have k^ t 



k 



at ■ 



i.e. the factorizable contributions satisfy our assumptions. 
We now consider the nonfactorizable contributions. If these pieces are inde- 
pendent of polarization (at least for the a = 0, || states), then in Eq. ()17|) we will 
have 

r = r\\=r , s = s\\ = s , p = P\\ = P , A = Ay = A . (23) 

This leads to k ct = k ut , so that our assumptions will be satisfied. Our assump- 
tions are therefore invalid only to the extent that the nonfactorizable pieces are 
transversity-dependent . 

In the heavy-quark limit with mb yC oo, there is only one universal form factor 
resulting from the spin symmetry of the theory. This implies that the (factorizable) 
Pg*' s i n Eq. (12 lj) are proportional for different polarization states. In other words, 
the various transversity amplitudes are related to one another. It is likely that 
these relations remain true in the presence of nonfactorizable corrections. This then 
implies Eq. (123)1 . We therefore expect deviations from Eq. (J23]l to be suppressed by 
0(1 /m c h ). In all, the net correction to our assumptions is 0{l/m c>b ) times the ratio 
of nonfactorizable and factorizable effects. We expect this to be small, so that the 
assumptions in Eq. (fTlj) are justified. 

Finally, we note that our assumptions can be tested. In the presence of nonfac- 
torizable effects of the form in Eq. (}2"3*|) . Eq. (|T^|) can be rewritten as 



A[B° d -> D*D*] 



G p 

V2 



X T Pp« , 



(24) 



where X T = Xvr + X 2 t, and 



X\T 

X 2 



V cb V* d a 2 + a\ + a[ — a% — a 



v ub v: d 



"10 



1 + re 



1 + se 



m D *f D *£% (D*\ C7 M (1 - 75)6 



B 



(25) 
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gamma (radians) 



Figure 1: The measured exclusion level, as a function of 7, from a fit to the 
vector-vector modes (D* + D*~ and D* + D*~). The exclusion level is defined in the 
text. From the fit, 7 is favored to lie in one of the ranges [0.31 — 1.50] (+0 or tt), 
[2.02 - 2.62] (+0 or tt), or [2.72 - 3.05] (+0 or tt) radians at 68% confidence level. 



As Xt is common to both a = 0, || states we find the relative phase between these 
transversity amplitudes is or tt. This prediction can be checked through an angular 
analysis of B — > D*D* decays. Note also that the assumptions in Eq. (fT4*|) are not 
in fact required in the original method of Ref. PJ. Eventually it will be possible to 
experimentally separate out the and the || components, making such assumptions 
unnecessary. 

With the assumptions of Eq. (fT4*j) . we can determine the value of 7 from the VV 
decays, using the method above. We use the measurements of the B — > £)(*)£)(*) 
and B — > D^D^ branching fractions, measurements of the B — > £)(*)£)(*) CP 
asymmetries, and the world-average values of sin 2/? (0.736 ± 0.049) [3 a and sin 2 6* c 
(0.0482 ± 0.0010) HE]. We take 2(3 to lie in the first quadrant. 

We use a toy Monte Carlo (MC) method to determine the confidence intervals 
for 7. We consider 500 values for 7, evenly spaced between and 2tt. For each 
value of 7 considered, we generate 25000 toy MC experiments, with inputs that 
span the range of the experimental errors of each quantity. For each experiment, we 
generate random values of each of the experimental inputs according to Gaussian 
distributions, with means and sigmas according to the measured central value and 
total errors on each experimental quantity. We make the assumption that the ratio 
Jd*/ Id* is equal to JdJ fo = 1-20 ± 0.06 ± 0.06 jS]. An additional theory error of 
10% is included to take into account the assumptions of Eq. (JHJ), as wen as smaller 
errors such as the neglect of the exchange diagram, subdominant SU(3)-breaking 
terms, etc. We then calculate the resulting values of Act, ^dir, ^Wir, and B, given 
the generated random values (based on the experimental values). Inputting the 
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Figure 2: A "pull distribution" check on the measured confidence levels for the 7 
fit. Values of 7, as well as the other theoretical parameters A ct , A ut , 5, and ft, are 
generated and used to produce values of the experimental inputs to the fit. The fit 
is then performed, and ((measured value of 7) — (generated value of 7)) /(measured 
uncertainty on 7) is plotted. The result is consistent with a Gaussian distribution 
with cr=l, implying that uncertainty on the measured value of 7 is accurately 
described by the confidence distribution. 



quantities a dir , a in di r , and B, along with ft and the value of 7 that is being considered, 
into Eq. (JBJ), we obtain a residual value for each experiment, equal to the difference 
of the left- and right-hand sides of the equation. One thus obtains an ensemble of 
residual values from the 25000 experiments. A likelihood, as a function of 7, can 
be obtained from x 2 {l)i where x 2 = (/V ") 2 ; i* 1 which /1 is the mean of the above 
ensemble of residual values and a is the usual square root of the variance. The value 
of xil) is then considered to represent a likelihood which is equal to that of a value 
X standard devations of a Gaussian distribution from the most likely value(s) of 7. 
We define the "exclusion level," as a function of the value of 7, as follows: the value 
of 7 is excluded from a range at a given C.L. if the exclusion level in that range of 
7 values is greater than the given C.L. 

Fig. ^ shows the resulting measured confidence as a function of 7. We see that 7 
is favored to lie in one of the ranges [0.31 — 1.50](+0 or 7r), [2.02 — 2.62](+0 or tt), or 
[2.72-3.05] (+0 or tt) radians at 68% C.L. This corresponds to [18.0° -85.7°] (+0° or 
180°), [116° - 150°](+0° or 180°), or [156° - 175°](+0° or 180°). Fig. H shows a 
check on the confidence distribution to ensure that it accurately describes the level 
of uncertainty on the measured value of 7. 

We now turn to the VP decays B° d (t) -> D* ± D T , B Q d -> D*+D~, and B° d -> 
D+D*~ . The advantage of the VP method is that no additional assumptions of the 
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gamma (radians) 



Figure 3: The measured exclusion level, as a function of 7, from a fit to the vector- 
pseudoscalar modes D* ± D =F , D* + D~, and D+D*~ . Unlike the vector-vector modes, 
with present data we do not obtain useful information on the likely regions of 7 from 
these modes alone. 



type described in Eq. flTIj) are needed. The disadvantage is that, as we will see, the 
data are such that no information on the most likely regions of 7 can be obtained 
from the VP modes. 

In order to implement the VP method, we proceed as follows. We first use the 
expressions for B, a dir , B and a^r in Eqs. © and (fTUj) to analytically solve for the 
theoretical unknowns A u t, A u t, 8, and 5', up to a four-fold ambiguity. We then have 
the two remaining equations for a^dir and a in( ar and two theoretical unknowns, 
7 and A. We now consider 200 values for each of 7 and A, each evenly spaced 
between and 2it. For each of the 200 x 200 possible combinations of values of 7 
and A, we generate 5000 toy MC experiments, with inputs that span the range of 
the experimental errors of each quantity. 

Similar to the toy Monte Carlo 7 determination for the VV mode, we generate 
random values of each of the experimental inputs according to Gaussian distribu- 
tions, with means and sigmas according to the measured central value and total 
errors on each experimental quantity. Making the assumption that fjj* / fu* is equal 
t° IdJ Id = 1-20 ± 0.06 ± 0.06 jS], we again obtain a confidence distribution as a 
function of 7. The result is shown in Fig. El As can been seen in this figure, present 
data on VP decays alone do not lead to useful constraints on 7. 

Finally, we can combine information from the VV and VP modes. The result 
is shown in Fig. |3J From the combined fit, we see that 7 is favored to lie in one of 
the ranges [0.34 - 1.41](+0 or tt), [2.09 - 2.57](+0 or vr), or [2.79 - 3.04](+0 or vr) 
radians at 68% confidence level. This corresponds to [19.4° — 80.6°] (+0° or 180°), 
[120° - 147°] (+0° or 180°), or [160° - 174°] (+0° or 180°). Comparing Figs. □ and 
0J we see that the favored ranges of 7 are slightly more constrained with the VP 
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Figure 4: The measured exclusion level, as a function of 7, from the combined 
information from vector-vector and vector-pseudoscalar modes. The combined in- 
formation implies that 7 is favored to lie in one of the ranges [0.34 — 1.41] (+0 or tt), 
[2.09 - 2.57] (+0 or vr), or [2.79 - 3.04] (+0 or tt) radians at 68% confidence level. 



data. Thus, although the VP data does not by itself constrain 7, its inclusion in the 
combined fit does have an effect. 

To summarize, we have presented the extraction of 7 using measurements of 
B° d (t) -> DW+DW- and B° d -> decays 0. We find that 7 is favored to 

lie in one of the ranges [19.4° - 80.6°] (+0° or 180°), [120° - 147°] (+0° or 180°), or 
[160° - 174°] (+0° or 180°) at 68% confidence level (the (+0° or 180°) represents an 
additional ambiguity for each range). The first of these ranges is that favored by 
fits to the Unitarity Triangle, assuming the standard model. The ranges come prin- 
cipally from data on vector-vector decays, although the vector-pseudoscalar decays 
do improve the constraints slightly. Note that, if we consider a larger confidence 
level, there are no constraints on 7 at present. However, this study demonstrates 
the feasibility of the method - with more data, we will be able to obtain strong 
constraints on 7. 
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